Electron backscattering from dynamical impurities in a Luttinger liquid 
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Electron backscattering in a Luttinger liquid with an impurity is investigated in the presence of 
zero point motion of the phonon lattice. The impurity can mean either a mass defect, an elastic 
defect or a pinning defect. The phonon spectrum is then affected by the presence of the defect, 
which enters in the renormalization group equations for backscattering. The RG equation becomes 
dependent on the energy cutoff for finite Debye frequency lod giving rise to finite energy effects. We 
compute the local density of states and show how the renormalization group flow is affected by a 
finite lud- 
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I. INTRODUCTION 

It is well known that electron-electron interactions give 
rise to Luttinger behavior in one-dimensional metallic 
systems. Standard studies assume the existence of an in- 
stantaneous interaction between electrons. Retardation 
effects need to be included if the interactions involve ex- 
citations, such as phonons, of energies much lower than 
the electronic bandwidth. An interesting class of materi- 
als where such excitations are expected to play a role are 
the carbon nanotubesP], where one dimensional features 
coexist with significant electron-electron and electron- 
phonon interactions. 

The effect of the coupling to acoustic phonons on 
a Luttinger liquid was studied in[2, Q. An interest- 
ing new feature is the existence of the Wentzel-Bardeen 
instability 0,1^ for sufficiently large couplings. The anal- 
ysis of the backscattering induced by an impurity on the 
transport properties of a Luttinger liquidQ has also at- 
tracted a great deal of interest. The extension of this 
problem to the case of a static impurity interacting with 
a coupled electron-phonon system is discussed inQ- 

In the present work, we generalize previous work to 
the case where the source of backscattering in a one di- 
mensional system has internal dynamics, and retardation 
effects have to be taken into account. We consider, in 
particular, the effect of a local modification of the elas- 
tic properties of a Luttinger liquid where the interaction 
between the electrons and the acoustic phonons is not 
negligible. We think that this study can be relevant for 
the carbon nanotubes. The electron-phonon coupling in 
nanotubes has been estimated to be significant (81], and it 
is possibly the origin of superconducting features at low 
temperatures [gl 1 10l| . For sufficiently small radii, the nan- 
otubes are expected to be close to the Wentzel-Bardeen 
instability The repulsive electron-electron interac- 
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FIG. 1: a) An impurity atom is inserted in the atomic chain, 
resulting in a mass defect, b) a structural deformation of the 
lattice results in an elastic defect, c) interaction of one atom 
with an underlying substrate results in a pinning defect. 



tion is also large [T^IT^IT3 |. Both impurities and phonons 
are expected to play an important role in the transport 
properties of nanotubes 15. 16, .ITL H^ . Note also that 
the Luttinger liquid characteristics of carbon nanotubes 
have been studied by measuring the changes induced by 
impurities and contacts on the transport properties _19y]. 

The next section discusses the model of the bulk sys- 
tem studied in this paper. The different types of defects 
considered here are described in section III. The method 
of calculation used, based on the Renormalization Group 
approach described in[^ is presented in section IV. A 
discussion and the main conclusions of the paper can be 
found in Section V. 



II. THE ELECTRON-PHONON CHAIN. 

Consider a 1-D chain (length L) of N atoms, mass to, 
each pair separated by a and elastic constant K between 
each pair. In the absence of inhomogeneities, the action 
describing the lattice is written in terms of the Fourier 
components of the lattice deformation field dk,uj- 
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where lud is the De bye fr equency. The sound velocity 
is c = auJu/2 — K /m, and /i = m/a stands for the 
hnear mass density. We wiU in the following take the 
continuum approximation 4sin^(fca/2)/a^ — > . 

The charge sector of the electronic degrees of freedom 
is described as a continuum Luttinger model (note that 
phonons and lattice deformations do not couple to the 
spin sector). Using a standard bosonized notation, we 
have: 
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Note the analogy between Eqs.Q and ^ in the limit of 
long wavelengths. One can identify fi with {uKp)~^ and 
the sound velocity c with the charge velocity u. Usually 
u is greater than c, since it is of the order oi vf in most 
cases, but this need not be the case in the presence of 
interactions. 

We describe the electron-phonon interaction by a de- 
formation potential: 
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(3) 



where the approximate sign holds for long wavelength 
excitations. The reduced coupling strength is defined as 
h = g^uKpl ^ and the reduced mass density as = 
fiuKp. 

The complete model can be then be expressed in terms 
of a diadic Green's function: 
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where the Green's function of this action then reads after 
inversion: 
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Here k = (ka)/2 is the dimensionless momentum, 
uj = uj/loo is the dimensionless frequency and v± = 



± y^Ac%^ + (c^ — u^)^) are the group veloc- 
ities of the polarons (hybridized electron and phonon 
modes) in the bulk. Note that the Wentzel-Bardeen in- 
stability comes about when v- becomes imaginary, i.e. 
for coupling strength b > u/c. This completes the de- 
scription of the homogeneous electron phonon system. 



III. THE DEFECTS. 



in the mass density at position xq is described by the 
action: 

UJ 

This results in a modified Green's function which is no 
longer translational invariant. In the absence of electron- 
phonon coupling, the phonon Green's function becomes: 
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with the notation Av^ = uKpAm/a, which ranges from 
Av^ = —v^ (the mass at position xq is zero) to infinity. 
The retarded, advanced, and causal Green's functions 
are obtained by adding a small imaginary part (with the 
corresponding sign) to and choosing the principal res- 
olution of the square root, in which the branch cut is in 
the negative real axis. 

In the absence of a mass defect, we obtain a simple 
phonon band for the local density of states. Introduc- 
ing a finite defect modifies the phonon spectrum. We 
can describe the changes induced by the defect by the 
transmission coefficient for phonons of energy w. The 
transmission coefficient at low energies is: 
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(8) 



Hence, lim^,^o 2^('^) — 1- The defect is tranparent for 
low energy phonons. This result is a consequence of the 
fact that the defect does not break the translational in- 
variance of the lattice. Hence, the system supports long 
wavelength phonons of low energy. The deviations from 
this low energy limit take place at energies w ~ m/ Am. 
For massive defects. Am 3> m, the transmission coeffi- 
cient tends to zero at finite, but low, energies. The lattice 
is effectively divided into two disconnected parts in this 
range of energies. 

The general case (6 7^ 0) requires some matrix inver- 
sion: 
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Analytic expressions can be obtained for all ranges of 
Am. Regardless of the presence of the defect, the electron 
band and the phonon bands hybridize, as the phonon 
LDOS shows a peak at frequency uud/c, similarly the 
electron spectrum shows a peak at lud- The fact that the 
perturbation enters in eq.© multiplied by is a sign 
of its irrelevance at low energies, as discussed earlier. 



A. Local change in the mass density. 



B. Local change in the elastic constants. 



We assume that the defect changes locally both the 
mass density and the lattice elastic constants. A change 



In the case of an elastic defect, one modifies the spring 
constant between two atoms of the chain xq and xi. The 
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action associated with this type of defect then reads 
1 
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Note that contrary to the case of a mass defect, here 
two phonon fields (at locations xq and xi) enter this ac- 
tion. In order to obtain the local density of states, we 
first derive the 3x3 local Green's function matrix Go as- 
sociated with fields 4>XQ,undxo,uj and d^i^uj in the absence 
of the perturbation: 



'^XQ .udxQ ,—UJ / 



^XQudxx —UJj 



I ^dxQ^uj4^XQ,~Uj) {dxQ,UjdxQ^—uj} {dxQ,Ujdxi^^ij} 

{dxi,Uj4^Xo,—Lu} {dxi,LudxQ^—Uj} {dxi ,Ludxi ^—iu} 



(11) 



which requires integration of the full Green's function of 
Eq. |SJl over momentum. Next, Dyson's equation is used 
to take the defect into account: 

G = Go + Go- AV(1- AVGo)-^ -Go (12) 

with A V the potential matrix associated with the action 
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In the absence of electron-phonon coupling, analytic ex- 
pressions for the perturbed phonon Green's function for 

dxQU) 5 



0122 are obtained: 
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(14) 

The phonon tansmission coefficient at low energies can 
be written as: 
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When the electron-phonon coupling is finite the imag- 
inary part of the local Green's function can be negative 
even if AK is within its allowed range. This result is a lo- 
cal counterpart of the Wentzel-Bardeen instability of the 
extended version. Similar instabilities have been studied 
in other strongly correlated svstems[20j. 

As in the case of a local change in the mass distribu- 
tion, a local modification of the elastic constant does not 
change the value of the transmission coefficient at low 
energy, as this perturbation does not alter the transla- 
tional invariance of the system. The range of energies 
for which the transmission coefficient is close to one is 

w ~ y'K/AK. 



C. Pinning of the lattice. 

We can also consider the case where the lattice is 
pinned by an external perturbation. Then, the action 
due to this defect is: 



^ / dtAKdl, 



ASk = -AV 

2 ^ ^ ^^dxQ—udxQU) 
Using the methods discussed earlier, we find: 
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In this case the perturbation is relevant at low energies, 
and the transmission coefficient goes to zero at low ener- 
gies. The lattice is effectively divided into two decoupled 
pieces. When the electron-phonon coupling is zero, the 
phonon transmission coefficient at low energies is: 



T{w) 



w 



w + i 



AK 
2K 



(18) 



Above a crossover energy, w > AK/K, the defect is 
transparent to the phonons, and the transmission coeffi- 
cient approaches one. 



IV. ELECTRON BACKSCATTERING. 
A. Flow equations. 

The defects considered in the previous section mod- 
ify the phonon and electron LDOS. Now we address 
how each type of defect affects electron backscattering. 
We consider here the limit of weak backscattering only. 
Then, without loss of generality, the action correspond- 
ing to a short range potential reads: 

ASe = yeJ2 f dr'f+.^^^^^rsxor = • ■ ■ = (19) 
rs 

= Se J dr cos (y2TT9pxor^ COS (y2TT9„xoT^ 

Following 7] , the flow of the electron backscattering 5e 
induced by the defect obeys the equation: 

A is the cutoff that scales as A = e^'Ap, and Q{xq,uj) is 
the Matsubara Green function, which results from per- 
forming a Wick's rotation on the action. After the Wick's 
rotation, using the a path integral formalism, the weight 
of the paths becomes e""^ instead of e*'^, where S is S — 
EkJ'l^,d]g-^[^:d]^ instead of 5 = EkJ'l^, d]G-^[cj,, d^ 



4 



Am 



Relevont 


Irrelevont 




1 

Unphysical 






FIG. 2: Backscattering flow diagram for a mass defect and 
two different values of tfie cutoff, A — ujd (left) and A — 
0.2ujD- This time the lower grey zone is unphysical, since 
it corresponds to a negative substitute mass. This time the 
defect resonance always remains in this region. 



8,(A) 







Relevant (k:<k„J 






— Irrelevant (K>K,,riJ 






Mass defect 



^^^^^^^^^^^ Pinning defect 

0,4 0.6 

A/co„ 



FIG. 3: Flow diagrams from a given initial backscattering 
term as function of energy. Top: defect which changes locally 
the mass with Am = 1. Bottom: defect which pins locally 
the lattice, with AK = 0.1. Different curves correspond to 
different k electron interaction parameters. Velocities were 
chosen to be m = 2c and b was set to a constant 2% below the 
Wentzel-Bardeen instability. 




FIG. 4: Backscattering flow diagram for an elastic defect and 
three different values of the cutoff, A = lud (left), A = 0.2ujd 
(right) and the case of very small cutoff A = O.OIljo (lower 
plot). The white and black regions denote relevant and ir- 
relevant flow of the backscattering respectively, while in the 
lower grey zone the system becomes unstable. The local insta- 
bility boundary (dashed region in the diagram) corresponds 
to AK/K < AKcrit/K = - (^^)^ In the lower plot we 
have emphasized the irrelevant backscattering sliver respon- 
sible for the effective suppression at flnite energies of elastic 
defect backscattering close to Ai^crit. 



(where G was chosen in the causal prescription for con- 
vergence). With these definitions it can be shown that 
Q{k,u!) = —G{k,iuj). This quantity is real, due to the 
fact that G{k, —oj) = G{k, uj). 

The main difference between the flow of the backscat- 
tering term in eq. H2U|) and the flow associated with elastic 
scattering in Luttinger liquids is the non trivial cutoff 
dependence in the right hand side in eq. . Retardation 
effects can also been induced by elastic defects in systems 
of electrons coupled to phononsQ due to the difference 
between the Debye energy, wd, and the electronic band- 
width. If the defect changes the elastic properties of the 
lattice, a different crossover at energies lower than lud is 
also possible, as discussed in section IIIB. 

At the lowest energies or temperatures, a local change 
in the clastic properties or the mass density does not af- 
fect the phonon transmission coefficient. At these scales, 
the defect plays no role, and the scaliiig of the backscat- 
tering term is given by the results inj^]- At higher en- 
ergies or temperatures, although they can be small com- 
pared to lod, the phonon transmission coefficient can be 
significantly reduced, more so the closer the system is 



to the bulk Wentzel-Bardeen instability. As the quasi- 
particles are made up of a combination of electron and 
phonon modes, this effect tends to enhance the electronic 
backscattering. 



B. Mass defects 

The sign of the flow of the backscattering term in 
ea. (|20|l allows us to divide the parameter space into re- 
gions where the flow is renormalized towards higher val- 
ues (relevant) or towards lower values (irrelevant) at a 
given energy. These regions are shown in Fig. j2] for 
a mass defect (see section IIIA). The separatrix between 
the two regions tends to be a straight line, independent of 
AM at low energies, in agreement with(7|, where a case 
equivalent to AM — was considered. At higher ener- 
gies, the region where the backscattering term appears to 
be relevant is enlarged. Hence, the flow of the backscat- 
tering term is not homogeneous in the regions where the 
right hand side of eq. H2U|I changes sign as function of 
energy. The flow of a given initial backscattering term is 
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FIG. 5: Backscattering flow diagram for a defect which elasti- 
caUy pins the lattice at a given site (due e.g. to the interaction 
with a substrate). 

shown in the upper part of Fig. The relevant flow 
is monotonous, although it deviates significantly from an 
exponential dependence on energy [3 near the critical 
line. The irrelevant region of the paremeter space, which 
is K > Krrit — "+"-('"++"-) jg nou-monotonous close to 
the WB instability, giving an enhanced backscattering of 
electrons when we are closed to the Kcrit boundary. 

C. Elastic defects 

The flow of the backscattering term is similar when 
the defect changes the elastic constants of the lattice (see 
section IIIB). The main difference is that a local Wentzel- 
Bardeen instability can take place even when the mod- 
ified elastic constants are positive, and the bulk is sta- 
ble. The boundary of the locally unstable region is given 

by ^ < ^ = - C-^r = -1 + {^y, and is 
marked by a dashed boundary in Fig. Near this 

Aifcrit hne we find a narrow sliver where the electronic 
backscattering term decreases at high energies, although 
it can eventually be relevant at low energies. This non 
monotonous behavior is the inverse of the one discussed 
for a mass defect. The asymptotic A ^ boundary for 
relevant-irrelevant behavior of backscattering on elastic 
defects is correctly derived in 0], and it is exactly the 
same as the Kcrit of the mass defect. However at finite 
temperature and k < Kent (repulsive electrons) backscat- 
tering from such defects can be strongly suppressed when 
we are close to Kcrit from below and to the local WB in- 
stability at Aif crit ■ The experimental relevance of this is 
greater than the bulk Wentzel-Bardeen instability, since 
the condition AK AK„it is much more easily achieved 
than the bulk instability, and can have measurable ef- 
fects on conductance at finite temperatures. Typical 
flow diagrams of backscattering amplitude for a given 




FIG. 6: Flow diagrams from a given initial backscattering 
term as function of energy. The defect changes locally the 
elastic constant. In the inset AK > while in the main 
plot AK was set to a negative values closer than 1% to the 
critical AK which is AK = —0.4375 for the parameters used. 
Velocities are u — 2c and b — 1.5 is a constant for all curves. 

initial value as function of energy and for different bulk 
parameters are given in Fig. |S] . One can clearly see how 
the non monotonous behavior discussed above is strongly 
enhanced near the critical line (main figure). 

D. Pinning of the lattice. 

A representative phase diagram for the case when the 
lattice is pinned by a defect which breaks translational 
invariance is shown in Fig. [S]. In this case, and at finite 
frequencies, the relevant region is reduced, giving rise to 
a non-monotonous renormalization of the backscattering 
amplitude that can be suppressed in a similar way as a 
mass defect could be enhanced at finite frequencies. This 
effect, as in the case of a substitute mass, is stronger the 
closer the system gets to the bulk instability. The flow of 
a given initial backscattering term is shown in the lower 
part of Fig. [3] for a system close to the Wentzel-Bardeen 
instability. Note moreover that the pinning defect case 
presents a peculiar scaling of the flow equations, which 
only depend on A/^ /A (see vertical axis of Fig. |3]). 

V. CONCLUSIONS. 

We have analyzed the role of lattice defects in Lut- 
tinger liquids with significant interactions between elec- 
trons and acoustical phonons. The quasiparticles of 
the system are combinations of electron-hole pairs and 
phonons. A lattice defect modifies locally the propaga- 
tion of the phonons, inducing a change in the electronic 
transport properties as well. 

A lattice defect which does not break the translational 
invariance of the lattice is irrelevant at low energies, as it 
is transparent to low frequency phonons. The borderline 
which separates the regions where electronic backscatter- 
ing is relevant or irrelevant is not affected by a defect of 
this type, and depends only on bulk parameters^ . There 
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is a range of energies, which can be much lower than 
the Debye frequency, ujd, where the defect scatters the 
phonons strongly. In this region of energies or temper- 
atures, the electronic backscattering can be significantly 
enhanced. Thus the flow of an electronic backscatter- 
ing term can be non monotonous, unlike the case of a 
Luttinger liquid with elastic scattering 6] , giving rise to 
measurable conductance dependence on temperature, for 
example. Even in a system with attractive interactions, 
a mass defect can block effectively the transport over a 
significant range of temperatures. On the other hand in 
a system with repulsive electronic interactions an elastic 
defect that is close enough to the local Wcntzcl-Bardcen 



instability can prove to be unexpectedly transparent to 
electrons. This local instability can be induced by a 
strong enough local softening of the elastic constant of 
the chain in the presence of electron- phonon interactions. 
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